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LUMINOUS SUPERNOVA-LIKE UV/OPTICAL/INFRARED TRANSIENTS ASSOCIATED WITH 
ULTRA-LONG GAMMA-RAY BURSTS FROM METAL-POOR BLUE SUPERGIANTS 
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ABSTRACT 

Metal-poor massive stars may typically end up their lives as blue supergiants (BSGs). Gamma-ray 
bursts (GRBs) from such progenitors could have ultra-long duration of relativistic jets. For example 
Population III (Pop III) GRBs at z ^ 10-20 might be observable as X-ray rich events with a typical 
duration of Tgo ~ 10'*(1 -I- z) sec. Recent GRB111209A at z = 0.677 has an ultra long duration of 
Tqq ~ 2.5 X lO'' sec so that it have been suggested that the progenitor might be a metal-poor BSGs 
in the local universe. Here, we suggest luminous UV/optical/infrared emissions associated with such 
a new class of GRB from metal poor BSGs. Before the jet head breaks out the progenitor envelope, 
the energy injected by the jet is stored in a hot-plasma cocoon, which finally emerges and expands 
as a baryon-loaded fireball. We show that the photospheric emissions from the cocoon fireball could 
be intrinsically very bright (Lpeak ~ io42-44 erg/sec) in UV/optical bands (epoak ^ 10 eV) with a 
typical duration of ^ 100 days in the rest frame. Such cocoon emissions from Pop HI GRB might be 
detectable in infrared bands at ~ years after Pop HI GRBs at up to z ^ 15 by up-coming facilities 
like JWST. We also suggest that GRB111209A might have been rebrightening in UV/optical bands 
up to an AB magnitude of < 26. The cocoon emissions from local metal-poor BSGs might have 
been already observed as luminous supernovae without GRB since they can be seen from the off-axis 
direction of the jet. 

Subject headings: stars: Population HI — gamma-ray burst: general — infrared: general 



1. INTRODUCTION 

Direct signals from Population HI (Pop HI) stars are 
needed to confirm their existence as well as to know the 
earliest history of star formation, galaxy evolution, and 
cosmic reionization at z > 10. However, to observe Pop 
HI stars in their stellar phase is extremely difficult due 
to the distance. Therefore the core-collapse phase may 
be more suitable to observe Pop HI stars. Especially 
possible supernovae and GRBs by Pop HI stars are good 
candidates. 

Now Pop HI stars are considered to be formed in a 
dark-matter halo of ~ IO^Mq at z > 20. The molecular 
hydrogen can cool the gas to ~ 200K, which determines 
the onset of t he Jeans in stability to form gas clum ps of ^ 
IO^Mq (e.g.. lBromm et a l. 1999: Nakamura fc Umemural 
I2002f ). If the most mass is accreted to form the star, 
the Pop HI star at the zero-age main sequence (ZAMS) 
would be as massive as Mzams ^ IO^-IO^Mq. If 
so, there have been discussed that Pop HI stars with 
ZAMS mass in a range of 140-260M(n would result 
in pa ir-instability supernovae (PISNe) tW ooslev et al.l 
120021) . Detailed numerical simulations have shown that 
the signals can be detect ed in infrared bands up to 
z ~ 30 using JWST (e.g. IWhalen et ai]l2012aD 0. As 
for even more massive stars, several core-collapse simu- 
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lations were performed, w hich would produc e black holes 
with - O(100)Mq (e.g. iFrver et alj [200lt iSuwa et al.l 
I2007allbh . 

However, such large ZAMS mass Pop HI stars seem to 
be formed only for spherically symmetric systems so that 
ZAMS mass of Pop HI stars may change when we take 
effects of rotation into account. Some cosmological sim- 
ulations indicate that the rotation of these gas clumps 
can natural l y spl i t them into sub c l umps of 5, lOOM in 
(iTurk et al.l[2009t IStacv et al.l [20101: IClark et al.l[20Tlh . 
Even if the seed gas clump is more massive, the protostel- 
lar UV radiation evaporates a significant fraction of the 
mass so that the Pop III stars are as small as Mzams 5 , 
(3O-9O)M0 (IMcKee fc Tanll2008l : lHosokawa et al.ll2011li . 
Population HI. 2 stars, which may b e more abundant, ca n 
also have smaller ZAMS mass re.g.. lYoshida et al.l[200l . 
Massive Pop III stars typically end up their lives as blue 
supergiants (BSGs) with massive h ydrogen envelope du e 
to the low opacity in the envelope ( Wooslev et al.|[2002|) . 

For such BSC Pop III stars, several authors 
have dis cussed the possibility of gamma-ray bursts 
(GRBs) JM eszaros fc Recs 2010; Komissarov & Barka 
^2UTd; 'Suwa fc loka 201 li iNagakura et all 120121 
iWooslev fc Hegeri ,2012.) . which we like to discuss 
from another point. One of the most promising sce- 
narios of long GRBs is collapsar scenario, in which 
relativistic jets burrowing through the stellar envelope 
are indispensable to bring relativistically moving ele- 
ments into optically thin region without dissipation. For 
massive progenitors with enough angular momentum, 
a stellar mass black hole (BH) with an accre tion disk 
would be first formed in th e core-collapse phase (jWooslevI 
119931 : iNaravan et al.|[2Q0ll ). and relativistic jets may be 
launched via the Blandford-Znajeck (BZ) like process 
(|Blandford fc Znaiekl [19771 IMcKinnev fc Gammiel[200l 
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iMeszaros fc Reed I2010D 0- If the jet can break out 
the stellar envelope before the central engine ceases, it 
may produce prompt gamma rays by some dissipation 
proc esses near and b eyond the photospheric radius (see 
e.g.. IMeszarosl l2012D . Also, one can expect radio to 
infrared afterglows via synchrotron emissions from the 
non-thermal electrons accelerated at the shocks in the 
relativistic shell going th rough the circumst e llar medium 
Jlpka & Mcszaros 200 51: llnoue et al.l 120071: iToma et al.] 
[201 It iNakauchi et al.ll2012f ).^ 

In terms of whether the jet can break out the enve- 
lope or not, the anticipated large pre-collapse mass of 
Pop III s tars w as thought to be problematic. Indeed, 
iMatzneJ ()2003D claimed that the progenitors of the ob- 
served long GRB would not be giant stars like red su- 
pergiants (RSGs) but compact stars like CO Wolf-Rayet 
(WR) stars. Nevertheless several authors have claimed 
that Pop III stars can produce GRBs irrespective of the 
pre-c o llapse mass ( Mcszaros & Rocs 2010; Suwa & loka 
iNagakura et al., .2012; .Wooslev fc Heger, ,2012i) . 
This is essentially because they would be BSGs in the 
pre-collapse phase (|Meszaros fc Reesll20?)l . The large 
mass gives a larger intrinsic energy budget, and the rela- 
tive compactness for the large mass helps the jet to reach 
the stellar edge before the engine ceases. As an interest- 
ing outcome, the duration of the prompt emission of Pop 
III or BSG GRBs would be much longer > 10^(1 -f- z) sec 
compared with the observed value (Suwa fc Iokal[20Tll : 
INakauchi eralll2012f ). INakauchi et al.l (|2012| ) argued the 
detectability of such ultra-long GRBs from BSGs and 
found that if the peak spectral energy and isotropic 
equivalent energy {Ep — Eiso) correlation holds the ob- 
served Ef"" ~ 100 keV with Ei^o ~ 10^^ erg and the 
duration Tgo = 6 x 10"* sec for Pop III BSG with mass 
40Mq at z = 9 which can be detected by an EXIST- 
type future instrument. These values are similar to 
the recent observation of ultra-long GRB111209A with 
Ei,o = 5.82 X 10^3 e rg, E"""^ = 310 keV and the duration 
of ~ 2.5 X 10"* sec (|Gendre et al.l[20Tl . Given that the 
observed rcdshift is relatively small, z = 0.677, the pro- 
genitor of this GRB might be a metal-poor BSG in the 
local universe. 

Here, we propose a new possible photon emission as- 
sociated with such ultra-long GRBs from Pop III stars, 
or more general BSG progenitors. We consider a cocoon 
i.e., a hot-plasma sheath of the jet, which is a neces- 
sary ingredient of the above collapsar model of GRB. 
Before the jet breaks out th e progenitor, the energy is 
stored in the cocoon (e.g., iLazzati fc Begelmanl 120051 : 
iBromberg et al.l[2011h , which can finally emerges and ex- 
pands as a baryon-loaded fireball. We evaluate photo- 
spheric emissions from such cocoon fireballs for various 
progenitors and discuss the detectability. We show that 
the emission is sensitive to the progenitor, so the cocoon 
fireball photospheric emission (hereafter CFPE) can be 
a diagnostic of Pop III and BSG GRBs. 

This paper is organized as follows. In SecISJ we intro- 

^ Another a lternative of the je t production mech anism is neu- 
trino process ijMacFadven fc "Woos lev 1999; Popham et al.l 119991 : 



duce GRB progenitor models which we use in the fol- 
lowing calculations. In SeclSl we model jet-cocoon for- 
mation and evolution before and after breaking out each 
progenitor. In SeclH we model CFPE, and evaluate the 
detectability. In SeclSJ we summarize our calculations. 

2. PROGENITOR MODEL 

We consider massive stars (Mzams > AOMq) with 
zero- {Z = 0), low- {Z = 10~'^Zq), and also solar metal- 
lici ty (Z = Zcd). We use p re-collapse stellar models given 
bv IWooslev et all (l200l . IHeger fc Wooslml (pfUOh and 
'Ohk ubo et al.l (|2009t ). which is listed in Table 1. The first 
column shows model names which include a letter rep- 
resenting metallicity (z-zero, 1-low, s-solar), the ZAMS 
mass, and the pre-collapse state (BSG or WR). The mass 
(M,) and the radius (i?,) in the pre-collapse phase are 
shown in the second and the third columns, respectively. 
One can see that only the progenitors with solar metallic- 
ity lose the dominant masses before collapsing to become 
WR stars. This is due to the strong stellar wind blowing 
out the hydrogen envelope, which may not be the case for 
zero- or low-metal stars due to the low opacity in the en- 
velope. Such mass losses also make the WR progenitors 
much more compact than the BSG progenitors. 

3. COCOON FORMATION AND EVOLUTION 

3.1. Before jet-cocoon breakouts 

In the case of massive progenitors with the pre-collapse 
mass > AOMq, a stellar mass BH would be f irst fo rmed 
just after the collapse is triggered (jHeger et alll2003[ ). We 
assume that the initial BH mass is 3Mq, and evaluate 
the time-dependent mass accretion onto the BH as 



M{t) = a 



dMr- 



dts 



t=ti!{r) 



(1) 



IZalamea fc Beloborodovl 12011k Istiwal 12013) . in which copious 
amount of neutrinos and anti-neutrinos emitted from the hyper- 
accreting disk annihilate in the vicinity of the axis and produce a 
fireball consisting of hot electron-positron pairs. 



Here = iirr'^ p^,{r')dr' is the mass coordinate for 
A'lr > 3Mq with being the density of the steller en- 
velope, tfi = {r^ /GMr)^^"^ is the free-fall time of a mass 
shell dMr with G being the gravitational constant, and 
a < 1 represents the deviation from the free -fall accretion 
13, which may be time-dependent (see e.g.. iKumar et all 
2008!). We use t as a time coordinate from the start of 
the accretion in the rest frame of the central engine. 

We assume that the luminosity of the relativistic jet is 
proportional to the mass accretion rate into the central 
BH; 

Lj(<)==r7M(^)c^ (2) 

where rj is the efhciency factor and c is the speed of 
light. E g. ([2t can be realized in t he BZ like pro- 
cess (e.g., iKomissarov fc Barkovl [20101 : iKawanaka et al.l 
[Ml . Note that rj may be also time-dependent. 

When the jet collides with the stellar envelope, two 
shocks are formed; a forward shock propagating in the 
stellar envelope, and a reverse shock in the jet. We 
call the region sandwiched by the two shocks as jet 
head. We consider jets with a Lor entz factor F; > 10 
(/3i c± 1). In the thin-sh eU limit (jSari fc Piran] 119951 : 
IMeszaros fc Waxmanll200l"h , which would be valid in our 



In ISuwa fcloE] 1 I20T1I) and INakauchi et al.l | |20T^ . the sup- 
pression factor of accretion rate, a, was absorbed in the definition 
of the jet efficiency 77; their ri corresponds to ar] in this paper. 
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TABLE 1 

Cocoon parameters at the jet breakout and the photon emissions at the photospheric radius. 



Progenitor 
model 


M, 

[M@] 


[em] 




[erg] 


Vc,bo 

[em'^] 


Mc,bo 
[Mq] 




£ 




^pcak 

[erg/see] 


^pcak 

[eV] 


^pcak 

[day] 


z40BSGt 


40 


1.5 X 10" 




1.8 X 10^2 


3.2 X 10^5 


0.92 


3.4 X lO'' 


0.011 




2.4 X 10*2 


13 


87 


zTOBSG* 


70 


1.4 X 10" 




3.9 X 10^2 


1.5 X 10^* 


0.88 


6.7 X 10^ 


0.025 




1.6 X 10"*^ 


11 


94 


z915BSGO 


915 


8.8 X 10" 




6.2 X 10^^ 


5.0 X IDS'" 


16 


2.2 X lO'' 


0.022 




2.1 X 10"*^ 


9.5 


350 


MOBSQt 


40 


4.4 X 10" 


fiducial 


1.8 X 10^2 


1.0 X 10=''' 


1.1 


3.9 X 10** 


8.8 X 10- 


-3 


4.8 X 10*2 


12 


91 


i75BSGt 


75 


8.6 X 10" 




4.5 X 10^2 


7.5 X IDS'" 


2.1 


1.9 X 10** 


0.012 




1.2 X 10"*=* 


11 


120 


s40WRt 


8.7 


2.3 X 10"' 




2.3 X 10^1 


6.5 X 1029 


0.11 


3.1 X 10" 


0.012 




6.8 X 10*" 


22 


40 


s75WRt 


6.3 


3.7 X 10"' 




1.1 X 10^1 


2.5 X lO^o 


0.075 


8.9 X 10" 


8.0 X 10" 


-3 


5.5 X 10*" 


21 


37 


/75BSG 






»7 X 10 


3.8 X 10^^ 


7.3 X IDS'" 


2.0 


1.9 X 10** 


0.10 




1.2 X 10** 


13 


68 


i75BSG 






r? X 1/10 


6.0 X 10^1 


1.2 X 10^8 


3.4 


1.9 X 10** 


9.7 X 10- 


-4 


1.2 X 10*2 


9.1 


220 


i75BSG 






ej X2 


5.3 X 10^2 


1.9 X 10^8 


5.3 


1.9 X 10** 


5.6 X 10- 


-3 


1.3 X 10*^ 


10 


140 


/75BSG 






X 1/2 


4.0 X 10^2 


3.4 X IDS'" 


0.96 


1.9 X 10** 


0.023 




1.1 X 10*^ 


11 


98 



The first column shows model names which include a letter representing metallicity (z-zero, 1-low, s-solar), the ZAMS mass, and the pre- 
collapse state (BSG or WR). The second and third columns show the mass and the radius of the progenitors, respectively. The 4th column 
shows the parameter of our theoretical model. Fiducial parameter means (a, ri, e) = (1, 6.2 X 10~*, 0.1). The 5th to 12th columns show the 
energy, the volume, the mass and the optical depth of the cocoon at the jet breakout, the baryon load parameter, the peak luminosity, the 
peak energy and the peak time of cocoon-fireball photospheric emission (CFPE), respectively. For details, see the text. 
The reference; ^ IWooslev et~all pOOl l. teeger fc Wooslev. 12010,), ^ .Ohkubo et al.l |(2509D 



cases Q, the two shocks approximately proceeds with a 
single Lorentz factor Fh- Then, the pressure balance 
at the jet head, (Lj/nefTj'^ri.^c) x (rjrh)2(l - AO^ = 
p,c^rh^/3h yields the velocity of the jet head as 



l + L(t)-i/2' 



where 



_ Lj(t-ri,/c) 



(3) 



(4) 



and 9j is the half opening angle of the jet. Here, the 
subscript "j" and "h" refer the jet and the jet head, re- 
spectively. The position of the jet head is obtained from 

rh = /o cl3hdt. 

As far as the jet head is non-relativistic, /3h ^ 
0.3, a dominant fraction of shocked plasma at the jet 
head will spread out sideways to form a cocoon (e.g., 
[Bromberg et al.ll20Tll ). The energy stored in the cocoon 
is evaluated as 



L;{t' -r^{t')/c)dt'. 



(5) 



Hereafter we simply assume the shape of the cocoon as 
circular cone with the height of ~ ''h(i)j s-i^d with the 
base radius of rc(<), where the subscript "c" refers to the 
cocoon. Then the volume of the cocoon can be evaluated 



as 



(6) 



The pressure balance at the cocoon-progenitor interface, 
Ec/Wc = p*/3c^c2, yields the expansion velocity of the 

The sufficient condition is L{t) <^ Fj*. 



cocoon as 



/3c (i) = 



1/2 



(7) 



The transverse size of the cocoon is given by rc = 
Jg cPcdt. Finally, one can evaluate the baryon mass 
loaded in the cocoon as 



4rh(t)= 



J 



(8) 



In principle, the cocoon is divided into an inner cocoon 
consisting of the shocked-jet matter and an outer cocoon 
of the shocked stellar matter. The two r egion can be sep- 
arate d by a contact discontinuity (e.g., IBromberg et al.l 
l20ll . In Eq. ©, we assume that the contact surface 
becomes unstable due to e.g., a Kelvin- Helmholtz type 
instability, and the outer and the inner cocoon are fully 
mixed. 

If the jet head arrives at the stellar surface before the 
central engine ceases, that is tho < ts{R*)/a, the jet 
succeeds to break out the stellar envelope, where tho is 
defined by rh(<bo) = R*- Hereafter the subscript "bo" 
stands for breakout. Note that the jet can successfully 
break out the stellar envelope for all the progenitor mod- 
els shown in Table 1. 

In our prescription, the three parameters a, ij, and 
9j determine whether the jet breaks out the envelope or 
not as well as parameters of the cocoon. In general, both 
a and rj depend on time, due to the details of the ac- 
cretion disk formation and the magnetic field structure 
around the BH. The jet-opening angle 6j also depends on 
time due to the comp ressio n by the cocoon pressu re (e.g., 
iLazzati fc Begelmanl[2005t IBromberg et al.ll20H"l ). In the 
following calculation, we use a = 1, 77 = 6.2 X 10-4, and 
0j =0.1 as fiducial values. By applying these to WR pro- 
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genitors, it is shown that the observ ed characteristics o f 
canonical GRBs can be reproduced (|Suwa fc Iokal[20TTI) . 

In Table m the internal energy (i?c,bo = £^c(ibo)), vol- 
ume (Vc^ho = K;(tbo)), and the baryon mass of the cocoon 
{Mc^ho = Vc{tho)) at the jet-cocoon breakout for vari- 
ous progenitors are shown. We also show the so called 
baryon-load parameter of the cocoon, 



Er 



Mc,boC2 ' 



(9) 



and the optical depth of the cocoon at the breakout, 



^c,bo 



m-nVc 



p t'c.bo 



(10) 



Here (Tt = 6.7 x 10 ^^cm^ is the Thomson cross section 
and 77ip is the proton mass. We should note that the 
cocoons are all non-relativistic^, i.e., e < 1, and also the 
cocoon fireballs are highly optically thick at the breakout, 
i.e., Tc^bo > 1- 

3.2. After jet- cocoon breakouts 

As shown below, the evolution of the cocoon after the 
jet-cocoon breaking out the progenitor can be character- 
ized by (i?*, i?c,bo, Afc,bo, V^.bo)- We do not consider an 
additional energy injection to the cocoon from the jet 
after the breakout for simplicity. 

The temperature of the cocoon fireball at the break- 
out can be determined from i?c,bo = a'7c,bo''K:,bo + 
-^^c.bo^B7c,bo/™p- In the case of cocoon fireball, the first 
term in the right-hand side is always dominant i.e., the 
fireball is radiation dominated, where the temperature 
evolution is described as Tc — Tc^bo (14/14, bo) "^^'^ with 



VaV'cbo 



1/4 



(11) 



Here a = 7.6 x lO^^^erg/cm^/K'' and fee = 1-4 x 
10~^^erg/K is the radiation constant and the Boltz- 
mann constant, respectively. Then, the radiation energy 
evolves as aTc^K « (14/K;,bo)~^^'^, and the gas energy 
evolves as (Mc_bo/"T.p)fcBTc oc (V;:/V;:,bo)~^^'^- Since the 
scaling is the same, if the radiation dominates initially, 
the situation will remain as far as the cocoon is optically 
thick. 

At first, the fireball would expand in an almost 
spherically symmetric manner, and the volume and 
the optical depth evolve as Vc = Vc^ho{R/R*)^, Tc — 
CTMc^hoR/mpVc = Tc^hoiR/ R*)'"^ , where R is the radius 
of the fireball. As long as the cocoon is optically thick, 
the kinetic energy Kc increases with the volume as dKc ~ 
p.dVc ~ aT,UVc/S = are,boVc^bo(i?/i?*)~'d(i?/i?,), 
which gives KdR) = -Ec,bo(l — (R/R*)^^). This means 
that once the cocoon fireball expands twice as large as 
the initial size, 

-Rsat ~ 2_R,, (12) 



^ In e.g., IPe'er et al.l (12003 '). I'ibma et aEI ||2009| '). and 
IStarline et"S] 1120 121) . the effects of relativistic cocoon fireballs have 
been discussed to explain the characteristics of observed GRBs. 
Cocoon fireballs can be relativistic if the mixing between inner 
and outer cocoon is suppressed. 



a considerable fraction of the internal energy will 
be delivered to the kinetic energy. The satu- 
ration velocity of the cocoon fireball is Wsat ~ 
(2£'c,bo/Mc,bo)^/^ = V2ec. The fluctuation of the 
velocity within the flreball 6v w (fcB?c/f^p)^^^ < 
(fcB?c,bo/27^p)^^^ is always much smaller than Wgat! 

Sv/vse^t < {Mc,hokBTc,holmpEc,hoY''^ < 1. Thus, the 
cocoon fireball can be approximated as a shell with 
a width A/i « i?sat beyond the saturation radius, 
where the volume and the temperature evolve as Vc oc 
(i?/i?sat)', Te cx (i^/i^,at)"'/^ respectively. The opti- 
cal depth still evolves as Tc = Tc^ho{R/ R*)~^ ■ Finally, 
the diffusion velocity of photons within the shell be- 
comes equal to the coasting velocity of the shell , that 
is, c/tc « Vsat- This occurs at the photospheric radius 
given by 

i?ph « 1.2i?, X (Tcbo'e)'/^ (13) 
In summary, cocoon fireballs evolve as 



c,bo 



(i?, < R < i?sat) 

2 (14) 



3 

(^) (^sat<i?<i?ph) 

r T A-i)" {R.<R<Rs..) 



R 



(16) 



4. COCOON-FIREBALL PHOTOSPHERIC EMISSIONS 

In this section, we evaluate the photon emissions from 
the cocoon fireball. Here we only consider thermal pho- 
tons and neglect other photon injection process, e.g. 
electron-synchrotron, bremsstrahlung emissions, which 
would give minor contributions in our cases. Also we ne- 
glect gamma ra ys as a decay product o f unstable nuclei 
like ^^Ni. In fact lTominaga et al.l (|2007[ ) showed that the 
abundance of ^^Ni synthesized by relativistic jets would 
be small. Moreover we neglect the effect of absorptions 
by nuclei included inside the fireball for simplicity. 

In the coasting phase, thermal photons within a width 
of « ARc/TcVsnt can escape the shell, since they have 
a larger diffusion velocity than the coasting velocity of 
the shell. For a fixed radius R, the radiation first comes 
from near the line of sight and later from the limb so 
that the photons from the fixed R are observed with a 
duration w R/vsm- Then, the mean bolometric luminos- 
ity can be evaluated as Lb ~ aTc"*^ x 47ri?^(Aflc/TcUsat) x 
{R/vsat)~^ (X (i?/i?sat)^^^ « i^/^, which becomes the 
maximum at the photospheric radius, R « Rph- The 
time scale of the emission is estimated to be ipcak ~ 
^ph/wsat as 



Vak = 0.84(i?,/c) X (TcboVe)'/*- 
The peak bolometric luminosity is described as L 

Ec ,bo X (Rsat/R*) ^(Rph/Rsat) ^^"^ X ^pcak 



(17) 



peak 



Lpeak = 0.84(i?c,boC/i?* 



X T 



-5/6 1/12 



(18) 
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Fig. 1. — Bolometric luminosities of CFPE from the progenitors 
in Table 1. We fix a = 1, = 6.2 X 10~^, and 9j = 0.1, and assume 
an exponential cut-off at t = tpcak- 



and the peak photon energy is £poak 

(i?sat/i?*)-^(i?ph/i?sat)-'/^ or 



3.92fcBT;,bo 



X 



£pcak = 2.8fcBTc,bo X (Tcbo'e)"'/^ 



(19) 



In Table[Il we show (ipoak, Epoak, ^pcak) for each progen- 
itor. 

Fig. [1] shows intrinsic bofometric luminosities of the 
cocoon- fireball photospheric emissions (CFPEs). We 
simply assume an exponential cutoff at at i = tpcak 
for each ca s e. Th is treatment would be justified since 
iPe'er et"al] (|2006f l numerically shows that the cutoff is 
steeper than oc for relativistic cocoon fireballs. The 
initial energy stored in the cocoon for z70BSG is about 
two times larger than that for z40BSG while the initial 
optical depth and the baryon load of zTOBSG are smaller 
than that of z40BSG (see Table [T|). Then, the cocoon of 
the zTOBSG case becomes transparent faster than that of 
the z40BSG case so that the former case is brighter. The 
initial optical depth for z915BSG is ^ 30 times larger 
than that for z70BSG while the baryon load e are simi- 
lar. Both the initial cocoon energy and the baryon mass 
of the z915BSG case are larger than those of the zTOBSG 
case. As a result, the evolution of the fireball looks sim- 
ilar, but the duration of the latter case becomes shorter. 

Comparing the progenitors with different metallicity, 
the CFPE from the metal-poor BSG progenitors with 
Z = 10^"^ Zq is similar to that of the Pop III progeni- 
tors, which can be understood from the fact that their 
pre-collapse stellar structures are also similar. In con- 
trast, the cocoon emission from the WR progenitors are 
much dimer than the other progenitors. This is essen- 
tially because that the timescale for the jet to break out 
the envelope is much shorter in the case of WR progeni- 
tors so that the energy stored in the cocoon is relatively 
small, about a tenth of the zBSG40 case. As a result, the 
peak luminosity is < 100 times dimmer. Therefore, the 
observation of the CFPE can be used as a diagnostic for 
successful Pop III GRBs and BSG GRBs. Moreover, one 
could constrain the progenitor mass from the observed 
parameters, (Lpoak, Epoak, ipoak), in principle. 

Next let us mention the dependence of the CFPE on 
the model parameters. Fig. [2] shows the cocoon emission 
from the /75BSG progenitor with different parameter sets 




Time after jet breakout [sec] 

Fig. 2. — Bolometric luminosities of CFPE from the i75BSG pro- 
genitor with different parameter set of (r], 6j) shown in Table 1. We 
fix o = 1, and assume an exponential cut-off at i = tpcak- 

of (rijO) with fixing a — 1 (see also Table 1). The solid 
line corresponds to our fiducial case, and the cases in 
which T] is 10 times larger (thick dotted line) or smaller 
(thin dotted line), and 6* is 2 times larger (thick dotted- 
dash line) or smaller (thin dotted-dash line) are shown. 
One can see that, by increasing r], the peak luminosities 
increase and the durations become shorter, and by in- 
creasing 9, the durations become longer, but the peak 
luminosities are almost the same. Using these informa- 
tions, we can probe (ry, 0), i.e., the efficiency and the half- 
opening angle of the jet inside the progenitor, from the 
observation of the cocoon emission once the progenitor 
is fixed, in principle. 

Now we investigate the spectral evolution of CFPEs 
in order to discuss the detectability using current and 
future facilities. The observed fiuxes in each frequency 
band can be described as 



T{t)Di^{zY 



(20) 



where tobs = ^(1 + z), Aobs = A(l -I- z), and Di,{z) is 
the luminosity distance of the source. We simply as- 
sume that the spectra are black-body ones, B\{T) — 
{2hc^/\^){exp{hc/XkBT) ~ 1)-^ for < R < Rpi, with 
h being the Planck constant. For a given wavelength, 
noting that the Rayleigh- Jeans law is applicable for op- 
tical and infrared radiation, ^\ob3(''obs) oc tobs^^^- In the 
engine rest frame, the fiux at wavelength A takes its max- 
imum at hc/X w 2.89fcBT'c(i) as long as it occurs before 
t = tpcak- This means that the observed flux takes its 
maximum at 



tobs, peak (Aobs) 

w min tpcak(l + -z), 



1.74t, 



(l + z)l/2 l^fcBTc^bo 



-3/2' 



(21) 



where tsat = -Rsat/^sat- 

Figs. 131 m and [5] show anticipated AB magnitudes of 
CFPE from Pop III stars in each band of optical to in- 
frared. We show the cases at z = 6 and 15 including the 
effect of Lyman-a absorption. Photons with observed 
wavelength Aobs ^ 0.85/xm and < 1.9/im are almost com- 
pletely absorbed in the case of z = 6 and 15, respec- 
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Fig. 3.— AB magnitudes of CFPE from a Pop III GRB with Mzams = 915Mq at z = 6 (left) and z = 15 (right). We fix a = 1, 
r] = 6.2 X lO"**, and 9j = 0.1. The horizontal line shows the anticipated lOcr-detection limit using JWST with an exposure time of lOOh. 
The effect of the Lyman-a absorption is included. 
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Fig. 4.— As in FigO but from a Pop III GRB with Mzams = 70M^ 
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Fig. 6. — As in Fig[2l but from a mctal-poor BSG GRB with 
4OM0 (top) and 75AfQ (bottom) at 2 = 0.7. 

tively, due to unionized interstellar medium. We neglect 
the e xtinction by intergalactic dusts (jWeinmann fc LilTvl 
l2005f ). The horizontal line shows the anticipated lOcr- 
detection limit using JWST NIRcam with an ex posure 
time of lOOh = 3.6 x 10^ sec0 (see also FZackrisso n et al.l 

unni). 

From Figs. |3] and HI one can see that JWST can 
detect the CFFE from Pop III GRBs with > TOAf© 
at up to z ^ 15. The anticipated duration would be 
longer than years. The event rate of such Fop III GRBs 
is still highly uncertain. Based on an optimistic as- 
sumptions, the all-sky event rate could be as high as 
< Q (IO^) yr~^ integrated over z > 6 (|de Souza et al.l 
[20T1I) . The FOV of JWST NIRcam with the shown sen- 
sitivity, 2.2 X 2.2 arcmin would not be wide enough for 
a blind search. Thus, a follow-up observation is needed, 
which is triggered by a wide-field X-ray telescope such as 
Lobst er and EXIST for th e prompt X ray and gamma 
rays (jNakauchi et al.ll2012D . or by a radio interferome- 
ters like ALM A, EVLA, LOFAR, and SKA for the ra- 
dio afterglow (l loka fc MeszarosI 120051 : llnoue et al.l 120071: 
[Toma ct al. 20n^^ 

The cocoon emissions from Fop III GRBs with < 
4OM0 would be relatively hard to detect at z > 6 (see 
Fig. [SJ. Nevertheless, such relatively low mass pro- 

^ http: / / www.stsci.edu/jwst / instruments / nircam/sensitivity/table 



genitors might exist even in a lower redshift, namely 
2 < 1. Fig. [6] shows anticipated AB magnitude from 
MOBSG and Z75BSG at z = 0.7. In these cases, an 
optical/infrared counterpart from months to a year af- 
ter the prompt emission would be expected. Interest- 
ingly, an observed ultra-long GRB111209A at z = 0.677, 
which exhibited very long duration of T90 ~ 2.5 x 10^ 
sec, might be origi nated from such a B SG progenitor in 
the local universe (jGendre et al.|[20T2l ). Based on our 
scenario, it might have been recently rebrightening in 
optical/infrared bands up to an AB magnitude of < 26, 
which could be detectable by current facilities. 

5. SUMMARY AND DISCUSSION 

We have investigated the possibility that hot-plasma 
cocoons associated with successful-GRB jets produce de- 
tectable signals. By using a simple model, we have cal- 
culated the formation and the evolution of cocoon before 
and after the jet breakouts, and the photospheric emis- 
sions from the cocoon for Fop III or metal-poor BSG, and 
also WR progenitors. We have found that the brightness 
can be highly progenitor-dependent, i.e., the bolomet- 
ric luminosity of the Fop III or BSG progenitors can 
be > 100 times larger than that of the WR progenitor. 
Thus, the cocoon-fireball photospheric emission (CFFE) 
can be a diagnostic of Fop III or BSG GRBs. We have 
shown that CFFE from Pop III GRBs even at z ~ 15 can 
be detectable in infrared bands as transients of years us- 
ing JWST. We have also considered the possibility that 
GRBs from metal-poor BSGs occur at relatively lower 
redshift, z < 1 [3- In these cases, our model predicts 
rebrightenings in optical/infrared bands from months to 
a year after the GRBs. Based on these results, we have 
proposed that an observed ultra- long GRB 111 209 A has 
been potentially followed by the CFFE recently up to an 
AB magnitude of < 26, which could be detectable. 

The local event rate of ultra-long GRB is estimated as 
- 9 X 10-3 Gpc-3yr-i (|Gendre et al.|[20T2l ) without tak- 
ing into account the beaming effect. Since CFFE can be 
seen from off-axis direction of the jet, the all-sky rate of 
the events can be as high as ~ 2(6'j/0.1)~^ Gpc~3yr~^, 
which is still much lower than a SN Ibc rate of ~ 
2 X 10^ Gpc-3yr"^ (|Madau et al.l 119981) . Nevertheless, 
they might have been already observed as a peculiar-type 
SNe given that the duration is long (< 100 days) and the 
total radiation energy is large (> 10^° erg). Especially, 
in the cases with a relatively energetic jet where rj is 10 
times larger than our fiducial value, total radiation en- 
ergy becomes as large as < lO^^erg (see Table 1), which 
can be comparable to that of the observed superlumi- 
nous supernovae (SLSNe). The e nergy source of SLSNe 
is still highly controversial ()Gal-Ya m 2012, and references 
therein). Given the local SLSN rate of ~ 10 Gpc~3yr~^, 
some of them might be originated from cocoon-fireballs 
associated with metal-poor BSG GRBs. 

In Fig. [71 we compare the light curve of cocoon emis- 
sions with that of various type SNe. We plot the bolo- 
metric magnitude of the cocoon emissions; the ^75BSG 
model with 10 times larger than our fiducial value for rj 
(thick solid line) and the s40WR model (thin solid line) . 
For the SNe, we plot the absolute magnitude in the ob- 



1" Actually, ri'ornatore et ai.l I I2007I ) showed that Pop III star 
formation could last up to z ~ 2. 
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Fig. 7. — Light curves of CFPE and SNe. The cocoon emission 
from the i75BSG progenitor with ?} = 6.2 X 10"^ and the s40WR 
progenitor are compared with jj = 6.2 X 10~^ are compared with 
SLSN-I (PTFOQcnd), SLSN-II (SN2006gy), SNIc associated with 
GRB (1998bw), and an average of the observed SNIbc. The cocoon 
emissions are plotted with the bolometric magnitude, and the SNe 
are with the R-band magnitude. 

served R band. We show a class I S LSN (PTFOQcnd; 
thick dotted hne; iQuimbv et all 120111) . a class II SLSN 
(SN2006gy; thin dotted line: iSmith et al]|2007[ ). a type 
Ic SNe associated wit h GRB or hyp ernova (SN1998bw; 
thick dotted-dash line: lGalama et al] il998). and an aver- 
age of the observe d type Ibc SNe (thin dotted-dash hne; 
iDrout et al.l[20Tll ). 

Prom Fig. [71 one can find that the CFPE from metal- 
poor BSG GRB can be similar to SLSNe at least in terms 
of the energetics and the timescale. In order to distin- 
guish between, or identify the cocoon emission and such 
a class of SNe, more detailed radiation transfer calcula- 
tions of the cocoon emission is required, which is beyond 
the scope of this paper. Fig. [7] also shows explicitly that 
the cocoon emission from the WR progenitor is much 
dimmer than the observed GRB SNe, or a typical type 
Ibc SNe. We can expect that the cocoon emissions as- 
sociated with WR GRBs are hard to detect since they 
would be hidden by possibly associated SNe or the host 
galaxy. 

Cocoon emissions are observationally characterized by 
three parameters (Lpoak, Cpoak, ipcak), from which one can 
obtain constraints between four physical parameters of 
the cocoon at the jet breakout (i?*, -Ec,bo, Afc,bo, V^,bo)- 
In our prescription, the latter three can be derived from 



the phenomenological parameters (a, ?7,0j). For a fixed 
parameter set of {a,r],9j), the observed characteristics 
highly depends on the progenitor mass, as shown in Fig. 
[TJ This means that one could constrain the mass of the 
Pop III or BSG progenitor by detecting the CFPE. On 
the other hand, if the progenitor mass is fixed, one can 
constrain (77, ^j) from the observed light curve (see Fig 12]). 
This means that one can potentially probe, e.g., the ac- 
cretion disk formation, the jet production and propaga- 
tion before the jet breakout, which otherwise cannot be 
seen only using electromagnetic signals of the prompt 
emissions. The observation of prompt gamma-ray and 
afterglow emissions from the relativistic jet could give a 
stronger constraint on (j],9j) and also a combined with 
the detection (or even non-detection) of cocoon emis- 
sions. 

Although our calculations are based on several simpli- 
fied assumptions, we can expect the results are qualita- 
tively correct. Our model parameters {a,r],9j), which 
are generally progenitor- and time-dependent, could be 
refined by numerical studies on the accretion-disk forma- 
tion and the jet propagation in the collapsing stars. In 
this paper, we assume that the inner and outer part of co- 
coon, i.e., shocked stellar matter and shocked jet matter, 
are fully mixed (Eq. ©). This estimate gives a maxi- 
mum baryon loading on the cocoons. If this is not the 
case, the cocoon emission would becom e bright er with 
a higher peak energy as discussed in Pe' er et all (|2006). 
The dynamics of cocoons just after the jet breakout, i.e., 
i?* < -R < i?sat, would be more complex than our model 
and can also be refined with the aid of numerical simu- 
lations. As for the coasting phase of the cocoon fireball, 
i.e., i?sat < -R < ^ph, more detailed radiation transfer 
calculations are needed to distinguish cocoon emissions 
from other competing sources e.g., various types of SNe 
and galaxies. 
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